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Abstract. In [WEI13b], the viscosity solution of the Hamilton- Jacobi equation was 
constructed by an "iterated minimax" procedure. Using Dafermos' front tracking 
method, we give another proof of this construction in the case of Hamilton-Jacobi 
equations in one space dimension. This allows us to get a better understanding in 
this case of the singularities of the viscosity solution. 



Consider the Cauchy problem for the Hamilton-Jacobi equations in one space variable 



There are in general no global classical C 1 solutions due to the crossing of charac- 
teristics. Different attempts to find proper "weak solutions" exist, such as Oleinik's 
and Kruzkov's entropy conditions, and explicit solutions constructed by Hopf formula 
[Hop65] for convex Hamiltonians or initial functions. 

The lack of uniqueness of weak solutions led M. G. Crandall, L. C. Evans, and P. 
L. Lions to introduce, in the 1980's, the notion of "viscosity solution" [Lio82, CEL84]. 
Viscosity solutions need not be differentiable anywhere, which makes their relationship 
with the classical crossing of characteristics unclear. However, they possess very general 
existence, uniqueness and stability properties and, in a large class of "good" cases, they 
coincide with the weak solutions introduced before. 

There is also a geometric method of constructing weak solutions for non-convex 
Hamiltonians, proposed by M. Chaperon, called "minmax solution", which, may not 
necessarily be the viscosity solution. In [WEI13b], the author has shown that, for 
general Hamilton-Jacobi equations, a limiting process of "iterated minmax" instead 
can lead to the viscosity solution. In one space variable, this fact can be explained 
more geometrically by investigating the wave front of the geometric solution. In this 
paper, we will introduce the method of "front tracking", which was first proposed 
by C. DAFERMOS ([Daf72]), to reveal the relation of iterated minmax with viscosity 
solution, and whence give an alternative perspective to the convergence of iterated 
minmax to viscosity solutions. As an application, in the last section, we will use 
the limiting process of iterated minmax to describe the propagation of singularities of 
viscosity solutions. 



1. Introduction 




d t u + H(d x u) = 0, 
u(0, x) = v(x), x e R 
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2. Preliminary on minmax solution 



We will briefly introduce the geometric frame work for Hamilton- Jacobi equations. 
We assume at first that H is C 2 and v is C 2 with bounded Lipschitz constant, denoted 
by Lip(v). The characteristics, or, the Hamiltonian flow of H, are given by 



where dv := {(x,dv(x)),x G M} is the 1-graph of v. Indeed, we can identify L with 
L = i(L) Cf(lx K) by the map 



A C 1 solution of the (HJ) equation, if it exists, is a function u(t, x) whose 1-graph 
{(t,x,d t u(t,x),d x u(t,x))} coincides with L, or equivalently, (x,d x u(t,x)) = tp t {dv) =: 
L t for all t. In general, due to the crossing of characteristics*, the geometric solution 
may be multi-valued, which prevents the existence of such a C l solution. 

The geometric solution can be generated by a C 2 family of functions St : K x M 2 — > M. 

(2.1) S t (x, x , y ) := v(x ) - tH(y ) + {x - x )y 

in the sense that, for all t, 



The function S(t, x, •) := St(x, •) : R x R x M 2 M is called a generating family of the 
geometric solution L. 

Note that L is contained in the subset U p := {\p\ < Lip(f)}. We can assume that 
H vanishes outside a neighbourhood of U p without changing L. If Q is the quadratic 
form in M 2 defined by Q(xq, yo) = —x Q y 0) then |'VS'(t, x, •) — VQ\ is uniformly bounded 
on compact subsets of (t, x). We say that S is a generating family quadratic at infinity. 
Indeed, for each compact subset of (t, x), S(t,x, ■) can be made exactly equal to the 
quadratic form Q outside a compact set by means of a fiberwise diffeomorphism. 

For any C 2 function / : X := IR fc — > R quadratic at infinity, such that f = Q outside 
a compact set, let f c := {rj\f(r]) < c} denote the sub-level set of /. Note that for 
c large enough, the homotopy types of f c , f~ c do not depend on c; hence we may 
denote them as /°° and f~°° . Suppose the quadratic form Q has Morse index A, taking 
the coefficient field Z2, the only non-zero homology groups are H\(f°°, f~°°; Z2) and 



tf n _ A (X\r-,X\/°°;Z 2 ). 

Let H and A be generators of H x (f°°, f-°°- Z 2 ) and H n _ x (X \ f-°°,X \ /°°;Z 2 ) 



V*(aJo, Vo) = (xo + tVH(y ),y ), 
and the geometric solution of the (HJ) equation is defined as 




i:Rx T*R -> T*(R x R), (t,x,p) ^ (t,x,-H(p),p) 



L t = ip\dv) = {(x,d x S t (x,xa,yo))\d ix0}yo) St{x,x ,y ) = 0} 



respectively. 



*It means precisely the crossing under the projection it : M x T*R 4ExI, (t, x.p) 1— > (t, x). 
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Definition 2.1. The minmax and maxmin of / are defined by 



inf max/ := inf max/(?7) 

[<j]=S r,e\cr\ 



min f := sup min 



where a is a relative cycle and |er| denotes its support. We call a a descending (resp. 
ascending) cycle if [a] = H (resp. [a] = A). 

The minmax and maxmin defined in such a way are equal, see [WEI11], and they 
are a critical value of /. 

Proposition 2.2. Assuming that v is a globally Lipschitz C 2 function and H is C 2 , 
the minmax 



is a weak solution of the (HJ) equation, i.e., it verifies the equation almost everywhere. 

We remark that since H is independent of t, the minmax R l s for any s, t depends 
only on the time difference t — s. The choice of the notation R\ instead of R l ~ s is 
simply for preference. 

The generating family quadratic at infinity defined by (2.1), and hence the minmax, 
can be extended to the Lipschitz framework, where v is globally Lipschitz and H locally 
Lipschitz, see [WEI13b]. In particular, when v is only Lipschitz, the initial 1-graph dv 
should be replaced by the enlarged pseudograph dv := {(x,p),p G dv(x)}, where d is 
Clarke's generalized derivative. 

Definition 2.3. Let / : IR fc — > R be a Lipschitz function, the Clarke's generalized 
derivative of / at x is defined by 



where co denotes the convex envelop. A point x is a critical point of / if G df(x). 

Example 2.4. For v(x) = \x\, dv(0) = [—1,1], the enlarged pseudograph dv is ob- 
tained by adding a vertical segment to the pseudograph dv := {(x, dv(x)),x G dom(v)}. 

The function S given in (2.1) is a generating family of the generalized geometric 
solution L = {J t {t} x cp^dv) in the sense that 

L t = (f t H (dv) := {(x,y )\y Q G dv(x ),x = x + tp } p G dH(y )} 
= {(x,d x S t (x,x o ,y o ))\0 G d( XOtyo )S t (x,x ,y )}, 

where (p* H : (x ,y ) i— > U p <=dH(y ){xo + tp,Uo) is the generalized Hamiltonian flow. We 
denote by C Lip (IR) the set of globally Lipschitz functions on R, by ||<9/|| the Lipschitz 
constant of / G C Lip (IR), and we denote \ f\n — max^g^ |/(^)| f° r an Y compact K C R. 
We rename by R t H v{x) the minmax function defined in(2.2), unless H is specified. 



(2.2) Rqv(x) := inf maxS^x, x , y ) = inf max 

(xo,yo) 




df(x) := co{limdf(x n ), x n -> x, x n G 



dom(d/)} 



4 QIAOLING WEI 

Proposition 2.5 ([WELL 3b]). Assume that H : R — > R locally Lipschitz and v G 
C Up (R), then, 

1) R l H v(x) is a critical value of the Lipschitz map (x ,y ) S t (x,x ,y ); 

2) T/ie minmax function R l H v(x) is Lipschitz and satisfies ||c?(i?^f)|| < ||<9t> ||; 

3) For any t\,ti > 0, 

|i?^t> (x) - i?#u(a;)| < \h - *2||-H"|{|p[<||&u[|}- 

4) Let H° and H l be two Hamiltonians, then 

\RhoV — Rtfivlco < \t\\H° — -f/' 1 |{| p |<||a 1 ,||}. 

5) If v ^ 1 G C Lip (IR) and K is a compact set in R ; then there is a compact set 
K = {x : \x\ < \x\k + T||<9if|{| p |<c}||} such that 

\tf H v° - tfjjV^K < l^-f 1 !^. 
where C = maxj ||<9t>j|| ; % — 0, 1. 

We recall the definition of viscosity solutions introduced by M.G.Crandall and P.L. Lions, 
which is defined, in general for first order partial differential equation, [CEL84]. 

Definition 2.6. A function u e C°((0, T) x R) is called a viscosity subsolution (resp. 
supers olution) of 

d t u + H{d x u) = 

when it has the following property: for every ip G C 1 ((0, T) x R) and every point (t, x) 
at which u — ip attains a local maximum (resp. minimum), one has 

d t i/) + H(d x ip) < 0, (resp. > 0). 

The function u is a viscosity solution if it is both a viscosity subsolution and superso- 
lution. 

For any convex Lipschitz funtion /, its convex conjugate /* is defined as 

/* : R — >> R = R U {oo}, f*(x) := sup(xy - f(y)). 

y 

Proposition 2.7 ([WEI13a]). If H G C 2 is convex, andv G C Lip (M), then the minmax 
is reduced to a minimum: 

Rov{x) = min (v(x ) - tH*( X ~ X ° )), t > 0. 

XO t 

This is one of the Hopf formulae which defines the viscosity solution of the (HJ) 
equation. There is another kind of Hopf formulae for convex inital functions v G 
C Up (R): 

(2.3) J*v{x) = {v* + tH)*{x) = m a x(xy - (v*(y ) + tH(y ))). 

yo 

It defines also the viscosity solution of the (HJ) equation. See [LR86]. 

In the following, we will denote by Jqv(x) the viscosity solution of the (HJ) equation 
with initial function v. 
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Lemma 2.8. Let v : R — >■ R be a convex Lipschitz function, and v* : R — > R be its 
convex conjugate, then 

y G dv(x) •<=>- v*(y) = xy — v(x). 

Proof. Since v is convex, its generalized derivative dv(x) is the usual subderivative, 

dv(x) = {y : Wx' , f(x') > f(x) + y(x' — x)} 
= {y : W, xy - f(x) > x'y - f(x')}. 

We conclude by the definition of the convex conjugate v*(y) = max^ x'y — v(x'). □ 

Proposition 2.9. We have the relation 

Rlv(x) < J l v{x). 

Proof. Recall that a generating family for the minmax is given by 

S t (x, xq, y ) = v(x ) + xy - tH(y ) - x y 

= (v(x ) - x y ) + xy -tH(y ). 

Let (x ,y ) be a point realizing the minmax: 

Rlv{x) = m£maxS t (x,x ,yo) = S t (x,x ,y ) 

Since the minmax is a critical value for the map (x ,y ) >->■ S t (x,x ,y ), we have 

G d Xo S t (x,x ,y ) y G dv(x ) <^> v*(y ) = x y - v(x ), 

hence 

Rlv(x) = xy - v*(y Q ) - tH(y ) 

and 

Rlvix) < m&x(xy - (v*{y ) + tH(y ))) = Jqv(x). 
yo 

□ 

In general, for convex (concave) initial functions and non-convex Hamiltonians, the 
minmax function and the viscosity solution may be different. The phenomena of rar- 
efaction serve as simple counterexamples. 

Example 2.10. Consider 

, . 1 6x/5, x < TT . . , 2 

V{x) = \-2x/3, x>0 ' H <P) = ~P +P +P- 

The initial function violates the viscosity condition, and there occurs rarefaction. For 
t > small and in a small neighborhood of x = 0, R t v(x) is not viscosity and Rqv(x) < 
Jqv(x). See Section 4 for a precise argument. 
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3. Front tracking and iterated minmax 

A notable feature of the viscosity solution is that, by its uniqueness, it possesses a 
semi-group property: 

Jqv(x) = J* o Jqv(x), 0<s<t. 

On the contrary, a minmax R l v(x) does not necessarily have this feature. One can 
always refers to rarefactions as counterexamples. The following proposition tells us 
that the semi-group property acounts exactly for the difference between the minmax 
and the viscosity solution. 

Proposition 3.1 ([WEI 13b]). The minmax Rqv(x) is the viscosity solution of the (HJ) 
equation if and only if it has the semi-group property: 

R t v(x) = R t s oR s Q v(x), 0<s<t. 

This motivates the construction of iterated minmax. Fixing a time interval [0,T], let 
Cn = {0 = t < ti < ■ ■ ■ < t n = T} be a subdivision of [0,T]. To each s G [0, T], we 
associate a number m(( n , s), depending on ( n : 

m(C„, s) := i, if U < s < t i+1 . 

For simplicity, fixing a subdivision, we may abbreviate m(( n , s) as m(n, s). 

Definition 3.2. The iterated minmax solution operator for the (HJ) equation with 
respect to a subdivision ( n is defined as follows: for < s' < s < T, 

K H,( n — K H ° ' ' ' ° K H 

When the Hamiltonian H is fixed , we may abbreviate our notation R s ^ := R l H s as 
Rl, and the iterated minmax as 

(3-4) Rl, n :=Rt, (n =Rt m(ns) o...R^\ 
which we call a n-step minmax, associated to a subdivision. 

Define the length of ( n by \( n \ := maxj \ti — t i+ i\. Suppose that {( n }n is a sequence 
of subdivisions of [0, T] such that \( n \ tends to zero as n goes to infinity, and let 
{i?Q n f(x)} n be the corresponding sequence of iterated minmax for an initial function 
v G C Lip (M). In [WEI13b], we have shown, for the general Hamilton- Jacobi equation, 
that any such sequence converges to the viscosity solution. 

For the (HJ) equation in one space variable, using the method of "front tracking", 
which was first proposed by C.DAFERMOS ([Daf72]), we can understand better the 
relation between iterated minmax and viscosity solution. 

Let us begin by considering the Riemann problem, with initial functions of the form 



v[x) 



P-X, x < 0, 
p + x, x > 0. 
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There are two possibilities: if p- < p + , then v is convex; if p_ > p + , v is concave. If 
P- < p + , denote by (H\[ p _ tP+ ^ the convex envelop of H on [p_,p + ], i.e. 

[p_,p+])"(p) ; — swp{h(p)\h < H, h convex on 

similarly, if p_ < p + , denote by (H\u +iP _])^ the concave envelop of H on [p + ,p_]. By 
the breakpoints of a piecewise linear function, we refer to the points where the function 
is not C 1 . 

By the Hopf formula (2.3), the viscosity solution of the Riemann problem of the (HJ) 
equation with piecewise linear H can be given explicitely. We take for example the 
convex case.: if (fl"|[p_ lP+ ]) w has m breakpoints in (p_,p + ), denoted by pi < ■ ■ ■ < p m , 
and p =p-, p m+ i = p+, set 

H(p i+1 ) - Hjpi) n ^ 

Si = , < i < m 

Pi+i - Pi 

Then 

{xp -tH(p ), x<ts , 
xp i+ i - tH(p i+1 ), x e [tsi,ts i+1 ], < i < m- 1, 
xp m+1 - tH(p m+1 ), x>ts m . 

In particular, Jqv(x) has m + 1 shock Xi{t) = s i^i < i < m. 

Remark 3.3. We remark that, by the use of the convex (resp. concave) envelop of H, it 
follows directly that, at each shock Xi{t) °f the viscosity solution Jqv(x), with the jump 
of derivatives Pi, Pi + i, the graph of H between p { and p i+ i lies above (resp. below) the 
segment joining (pi,H(pi)) and (pi+i, H(p i+ i)). This is the so-called Oleinik's entropy 
condition for viscosity solutions, see Lemma 4.1 in the next section. 

Now we are willing to investigate the minmax solution under the same hypotheses. 
We first give a profile for the wave front. The wave front of the geometric solution at 
time t is given by 

T l := {(x, S t (x,x ,y ))\yo G dv(x ),x G x + tdH(y )}. 

For any subset 4cK, define 

P A := FlixoeA} = {(x,S t (x,z )) G J=\x$ G A}. 

We claim that the wave front T 1 for ^(dv) with v and H piecewise linear (with 
finite pieces) is formed by pieces of straight line segments. Indeed, 

1) T\ := F\ X0>Q } and T^L := ^F\ X(I<Q ^ are two lines with slope p + and p- respectively. 
Take the case x < for example, one has y = v'(xq) = p_, and 

F{xo} = i z ( x 0i V) = ( x o + ty, v(x ) + t(yp- - H(p-))) : y G dH(p_)}. 

Then, for any x , x' < 0, y, y' G dH(p_), the chord connecting z(x ,y) and z(x' ,y') 
is of slope p_ . 

2) Without loss of generality, we assume that p_ = p' < ■ ■ ■ < p' k = p + (orp + = p' < 
" ' ' < Pk = P-) are the breakpoints of H between p_ and p + . Then = J rt \{ Xo =o} 
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consists of k line segments with slope p\ which correspond to the breakpoints p\. This 
can be seen from the formula 

^{0} = i z (y^p) = (ty,t(yp- h (p))) '■ v e \p+,p-],ye OH(p)}. 

If p + , p- are not breakpoints of H, then J~^ y loses two line segments of slopes p + and 
p_, but the whole wave front J 7 ' does not change in the presence of the two segments 
J 7 ! and T\. 

Example 3.4. Let p- < p + , which relates to a convex Riemann initial function v. 
Figure 1 gives H and its corresponding wave front. The minmax is obtained by choosing 
the max. 




Figure 1. wave front 



Lemma 3.5. For the Riemann problem with piecewise linear Hamiltonian, suppose 
that the graph of H between p- and p + lies above (resp. below) the segment joining 
ip-,HipJ)) and (p + , H(p + )) . Assuming p- < p + (resp. p_ >p + ), then J-t \ lies below 
(resp. above) the graph of the viscosity solution Jqv{x) in the wave front. 

Proof. We will give the proof in the case where p- < p+, while the other case is 
similar. For convenience, we may assume that p± are not breakpoints of H. The 
viscosity solution Jqv(x) has a shock x(t) at which J 7 ! and T\ intersect, 

= H{p-)-H{p + ) t = + = + 

V- ~P+ 

for some < and x + (t) > 0. We will show that, at x — x(t), 3~\o} lies below 

(x(t)iP— x-(t)+t(p— H'(p-)— H (p_))). The points in J 7 ^! are given by (ty,t(yp—H(p))) 
with y G dH{p) and p breakpoint of H in (p_,p + ). Let ty = x(t) — x -(t) + tH'(p_), 

t{yp - H{p)) - (p_z_(t) + t(p_ff'(p_) - H{p J))) 

= t(yp - Hip)) - t(yp-Hip-)) 

Hip) -Hip J). 

= tip-p.){y--^ ^ <0, 

p — p- 
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where the inequality comes from y = H ^ p p and the hypothesis on the graph of 

H. Hence J^ y lies below the graph of Jqv(x). In particular, they can intersect only at 

x(t). " " ' □ 

Definition 3.6. We say that a G K. admits a descending (resp. ascending) cycle if there 
is a descending (resp. ascending) cycle a along which a is the maximum (resp. minimum) 
of the generating function S. 

Lemma 3.7. //aeR admits at the same time a descending cycle and an ascending 
cycle, then a is both the minmax and maxmin value. 

Proof. By definition, it is easy to see that inf max S < a < sup min S. The inverse 
inequality follows from the fact that a descending cycle and an ascending cycle must 
intersect. □ 

Proposition 3.8. For the Riemann problem with piecewise linear Hamiltonian, the 
minmax solution Rqv(x) coincides with the viscosity solution. 

Proof. We first remark that for an arbitrary initial function, the minmax and the 
viscosity solution may differ immediately. 

Consider the Riemann problem with initial value v{x) = p-X, x < 0, v(x) = 
p + x, x > 0, with p_ < p + . It is sufficient to prove that Rqv(x) is piecewise lin- 
ear and has m + 1 shocks x«(0 — tsi, < % < m, where Sj = with 
P- = po < ■ ■ ■ < p m +i = P+ the breakpoints of the convex envelope 

For a fixed time t, let Vi, < i < m, be the intersection point of the line segments 
corresponding to Pi and p i+ i in the wave front J 7 ' 

Vi = (tSi,t(siPi - H{pi))) = (tsi,t(sip i+1 - H(p i+1 ))). 

We want to show that the minmax Rqv(x) is obtained from J 7 by selecting the 
segments ViV i+ i, < i < m — 1, 



ViV i+1 = {(ts,t(sp i+1 - H(p i+1 ))) : s between s { and s i+1 }. 

Fix any i, < i < m — 1, note c| := t(sp i+ i — H(p i+1 )). We claim that (ts, cf) admits 
a descending simplex, i.e. there exists a descending simplex a? such that 

max S t (ts, z ) = C*. 

Write explicitly 

S t {ts, z Q ) - c\ = p + Xq - tH(y ) + (ts - x )y - t(sp i+1 - H(p i+1 )) 

= Mp+ - Vo) + tA i+ljS (y ), if x > 0. 
S t (ts, z ) - c\ = x (p- - y ) +tA i+1;S (y ), if x < 0, 



where 

A+i,s(yo) 



\vo ~ Pi + i) (s ~ ^gg^ ) , for yo^Pm, 
0, otherwise. 
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By the convexity of (H |[ p _ iP+ ]) w ', we get 

A i+ i >s {y ) < 0, if y E \p-,p+), s between and s i+1 . 

Note that the minmax and maxmin depend on the Hamiltonian H only on the compact 
set \p-,p + ]. For each fixed i, we can modify H outside [p_,p+] by requiring 

H(y ) - H(p i+1 ) = f HiP ^r\ Vo<P-, 

so that 

A i+ i,s{y ) < 0, Vyo e R. 
We construct a descending simplex erf = (£0(2/0X2/0) as follows: 



£0(2/0) 



+ %o), Vo e [Pi+i,oo) 



t^z+i,s(yo) 
2/0 -p- 

^±^-%o), 2/o G (-oo,p i+lj 



where 9 : R — >■ R is a non negative continuous function, with 



«(yo) 



0, yo G [p_,p+], 
ll/ol, |l/o| > C, 



with C > a large constant. One sees that a\ is a descending simplex since the first 
term in Xo(yo) is bounded, thus 



, !i m 1 2/o I 1 ((2Co(l/b),3/o) - (2/o,2/o)) = 0. 



Furthermore 



6"(2/o) : = (zo (2/0), 2/o)) - c| = 0, if y G [p_,p+] 

For y e 00), the first term in £0(2/0) is negative, and one can verify that &f (2/0) is 
negative both in the cases where £0(2/0) < and £0(2/0) > 0. Similarly, we can show 
that bf(yo) < for y e (— 00, Therefore c| admits erf as a descending simplex. 

On the other hand, replacing 6 by —6 will give us an ascending simplex for cf, i.e. 
bj,(yo) > 0, with equality for y G [p_, Hence is at the same time a minmax and 
maxmin value by Lemma 3.7. 

For the case where p- > p+, that is, when v is concave, we should take the concave 
envelop (H\^ p+tP+ ])'^ and the proof is similar. □ 

Now suppose that v is piecewise linear continuous (with finite pieces). Thanks to its 
local nature and semi-group property, one can construct the corresponding viscosity 
solution by viewing v as a combination of Riemann initial data and chasing the inter- 
actions between the shocks of each sub Riemann problem. This is the so-called front 
tracking method, which consists in the following inductive procedure: every time there 
are collisions between the shocks, we restart by considering the resulting state as a new 
initial function and propagate until the next time of collision. In the presence of finite 
propagating speed of characteristics, the number of shocks of the consequent solution 
will never blow up. 
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Proposition 3.9 ([HR11], Lemma 2.6). Assume that v(x) is a piecewise linear contin- 
uous function with a finite number of discontinuities of dv(x), and H is a piecewise 
linear continuous function with a finite number of breakpoints in {\p\ < \\dv\\}. Then 
the viscosity solution of the (HJ) equation u(t,x) is piecewise linear in x for each t, 
and d x u(t,x) takes values in the finite set 

{dv(x),x G dom(dv)} [^{breakpoints of H}. 

Furthermore, the number of shocks of u is finite, bounded by a constant uniformly for 
all t, and the number of collisions of shocks is finite. 

Lemma 3.10. Suppose that t > 0, and O C IR is an open subset. If for any e > 0, 

there exists s G (0, e) such that Rl o R^v(x) is a viscosity solution of the (HJ) equation 
for (t,x) G (s,t) x O, then Rqv(x) is also a viscosity solution for (t,x) G (0,t) x O. 

Proof. We will apply the standard argument for the stability of viscosity solutions. Let 
s n > be a sequence of decreasing numbers such that s n — > 0, and note 

u n (r,x) := R T Sn oR^v(x), (t,x) G Sn := (s n ,t) x O 

and u(r, x) := Rqv(x). 

Given any y := (r, x) G O , let ip G C 1 (O f / 2 ) such that u — ip attains a local 
maximum at (f,x). Take ip' G C 1 {O f /2) such that < tp' < 1 if (r, x) —: y ^ y and 
ip'(y) = 1. Then u — (%p — ip') attains a strict local maximum at y. For n large enough, 
u n is well-defined in Of/2, and by Proposition 2.5, u n — > u on O f /2, thus there exists 
y n G Of/2 such that u n — — ip') attains a local maximum at y n and y n — > y. By 
assumption, u n are viscosity solutions, hence 

dtty - il>')(y n ) + H(d x (^ - *l)'){y n )) < 0, 

and we conclude that 

d t iP(y)+H(d x i/;(y))<0, 

since dip'(y) = (d t ip', d x ijj')(y) = 0. Thus we have proved that u is a viscosity subso- 
lution. Similarly, we can prove that it is a viscosity supersolution, hence a viscosity 
solution. □ 

Proposition 3.11. Assume that H and v satisfy the hypotheses in Proposition 3.9. 
If T\ > is such that for t G (0,Ti), there is no collision of shocks of the viscosity 
solution of the (HJ) problem, then the minmax solution R t Q v(x) coincides with the 
viscosity solution for t G (0, Ti] . 

Proof. Let {xj}i<j< ni be the set discontinuities of dv, and Xjif), 1 ^ 3 ^ n 2 be the 
shocks of the viscosity solution Jqv(x) for t G (0, Ti), where Xjif) < Xj+i(t)- Denote 
by J-^iv) the wave front at time t originated from v. Let e > be small enough 
such that the sub Riemann problems of initial function v are independent in the sense 
that those J~f x .\(v) emitted from dv(xi) in the wave front have no intersections with 
each other for t G (0, e). Then for t G (0, e), by the local property of minmax solution, 
R t Q v(x) is chosen simply by combining every sub minmax problem, hence it is piecewise 
linear, and by the local property of viscosity solution, it is viscosity. For any s G (0, e), 
writing v s (x) := Rqv(x), we claim that u Sjt (x) := R^oRqV^x) = Rlv s (x) is the viscosity 
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solution for t G (s,Ti). For t — s < €\ small enough, when there is no intersections 
between the T\ 



t-s 



{xj{s)Y 



, 1 < j < 712, we can apply the same argument as before : the 
mmmax u Si t(x), which is obtained by taking the minmax of each sub Riemann problem 
with the singularity Xj( s )i is the viscosity solution ; In particular, by Remark 3.3 and 
Lemma 3.5, every J 7 ! 



t-s 



is no collision of shocks, J 7 ! 



t-s 



lies above or below the graph of u s j(x). As long as there 

can intersect only strictly above 



and ^ { Ws)}^ 



or below the graph of u Sjt (x). As the minmax u s (t, x) is continuous in t, it cannot 
instantaneously jump to the intersection, hence it preserves the shocks Xj{t)> thus is 
the viscosity solution. 

This proves that, for e > small enough, and any s G (0, e), R\ o Rqv(x), t G {s,T\) 
is the viscosity solution. By Lemma 3.10, as s — > 0, R\ o Rf j v(x) — > Rqv(x) which is 
the viscosity solution for t G (0, T\]. 

□ 



Corollary 3.12. Assume that H andv satisfy the hypotheses in Proposition 3.9. Then 
there exist k > and Tj, 1 < j < k with < T\ < T 2 < ■ ■ ■ < such that, 



(3.5) 



Jlv(x) = Ef r .o---oI&v(x), Tj<t<T j+1 , 0<j<k. 



where, by convention, we write T = and T^ +1 = +oo. Furthermore, for any sequence 
of subdivisions {( n }n such that \( n \ — > 0, the sequence of iterated minmax {Rq ^ n v(x)} n 
converges to the viscosity solution 

Proof. Let Tj, 1 < j < k be the times of collisions of shocks of Jqv(x). By Proposition 
3.9, k is finite and we obtain (3.5) by Proposition 3.11. For the second part, given a 
T > T fc , let ( = {1 < ti < t 2 < ■ ■ ■ < t n = T} be a subdivision of [0, T}; we will show 
that 



(3.6) 



\Rl,Cn V - J ^lco < 2A;|#| {b |<| W} |C|, t G [0,T]. 



Denote by £' the subdivision obtained by adding to ( the Tj, 1 < j < k, as new 
dividing points. Then by the semi-group property of the viscosity solution, we have 

: maxi<j< n {/f:j, ti < Tj} and notice that for \(\ small, 



R^vix 



Jqv(x). Let Tj 



we have Tj < Tj +l . For t G [0, Ti], R t {)C y{x) = Jqv{x). Writing Vj(x) := i? ^f(a;) and 
v'j(x) := R ^,v(x), then for t G [T j+1 , T j+2 ], 





Kc v - 




1 R* 


< 


1 f?* 

'^- +1 ,c 


+ 


1 f?* 


< 


\Rf\ 


< 


\Vj - V'j 



o R T J+1 Vj - R\ 



o o R 



o rZ s+1 



o R T > +1 



Vj - R\ 



o i?^- 1 ° n 

T i+i 



+ Y|C0 



+ Y|co 



'j'+l 

^"+1 



j+1 v' 



T j+l\\ H \{\p\<\\dv\\}- 
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Hence, by induction on j, we can obtain (3.6). It follows that, if {Cn}n is a sequence 
of subdivisions of [0,T], such that \( n \ —> 0, we have 




Theorem 3.13. Suppose v : R — > R is globally Lipschitz and H : R — >■ R locally 
Lipschitz, then for any sequence of subdivisions {( n }n of[0,T] such that \( n \ — > 0, the 
sequence of iterated minmax {R°^ n v(x)} n converges uniformly on compact subsets to 
the viscosity solution of the {HJ) equation 

f d t u + H(d x u) = 0, te(0,T} 
1 u(0,x) = v(x), x G R 

Proof. For any compact subset if C 1, fixed once for all, we can modify if and v 
so that ii is linear outside K := < ||<9i>||} and v is linear outside K = {\x \ < 
\x\k +T\\d(H\jc)\\}. For k G N, take the piecewise linear approximations of t> fe and H k 
for f and ii such that 

\v - v k \ c o < \\dv\\/k, \H - H k \ c o < \\dH\\/k. 

Since H and v are linear outside compact sets, we can assume that both v k and H k 
have finite pieces. 

Given any sequence of subdivisions {( n } of [0, T] such that \( n \ —> 0, let P^n n v [x] be 
the corresponding iterated minmax. It follows from Propostion 2.5 that the sequence 
u n (t,x) := RffnV^x) is equi-Lipschitz and uniformly bounded for (t,x) G [0,T] x K. 
Hence by the Arzela-Ascoli Theorem, {u n } n has some convergent subsequence. We 
denote by R l H v(x) its limit. We have 

\tf Htn v - # H k tn v k \ K <\v-v% + T\H - H k \ K . 

By Corollary 3.12, for each k, the sequence {R t H k , n v k\n converges as n — > oo; denoting 
its limit by B? Hh Vk(x), it is a viscosity solution of the equation d t u + H k (d x u) = 0. We 
have 

\& H v - Rh^Ik <\v- v k \ k + T\H - H k \ K ^0, k oo. 

We get from the stability of viscosity solutions (ref.Theorem 1.4 [CEL84]), that Ftf H v{x) 
is a viscosity solution for d t u + H(d x u) = 0. By the uniqueness of the viscosity so- 
lution for the Cauchy problem, the limit R l H v(x) is independent of the convergent 
subsequence. Thus we conclude that the sequence of iterated minmax {R^ n v(x)} n 
converges to the viscosity solution. □ 



4. Singularities of viscosity solution 

In this section, we will see how their interpretation as a limit of minmax help to 
describe the propogation of singularities of viscosity solutions. 

The viscosity solution of the (HJ) equation is characterized equivalently by Oleinik's 
entropy condition: 
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Assume that O C (0, oo) x R is an open set, u G C(0), and that there is a C 1 curve 
X dividing O into two open subsets + and 0~, O = + U x U 

Lemma 4.1 (Entropy condition.). Lei w G C(0) and u = u + in + U X; u — u ~ i n 

0~ Ux, where G C 1 (0 ± Ux). T/ien u zs £/ie viscosity solution of the equation (HJ) 
in O if and only if: 

(1) u are classical solutions in O respectively, 

(2) The graph of H lies below (resp. above) the line segment joining the points 
(Pt, H (Pt)) a nd(pt,H(pt)) ifpj < pi (resp.pt < pf), wherepf := 9 a; M ± (t, 

In particular, we say that x(t) strictly verifies the entropy condition if it verifies the 
entropy condition and the line segment joining (p^ , H{p^)) and (pf,H(pf.)) is not 
tangent to the graph of H. 

Remark 4.2. Equivalently, the entropy condition can be described, for example, when 

Pt < PT, as 

j/W-gfa-) > gfrft-gfa-) 

w - Pt Pt - 

for any iy between p^ 1 " and The "=" between brackets refers to the Rankine-Hugoniot 
condition which is a necessary condition for u to be a weak solution in the sense of 
distribution. 

In the following, we suppose that the Hamiltonian H is smooth, with generic critical 
and inflection points, and that the initial function v is globally Lipschitzian. 
Recall that the wave front of the {HJ) equation at time t is 

JF = P(v) := {(x,S t {x;x ,yo))\yo G dv(x ), x = x + tH'(y )}, 

where S t (x, xq, yo) = v(xq) + xyo — tH(y ) — x 2/o is the generating family. It is a curve 
in the plane M 2 having cusps and self-intersection points. Adopting the notions in 
[EH87], a wave front is said typical if it has no triple self-intersection points and the 
sets C and D of cusp and double self-intersection points respectively are discrete and 
disjoint. The branches of the front J 7 are the irreducible components of J 7 \ C. Every 
branch is the graph of a partially defined C 1 curve. The union in space-time of the 
fronts is called the big front, denoted by J 7 , 

f = T(v) := {(t,x,St(x;x ,y ))\y Q G dv(x ),x = x + tH'(y )}. 

The projection of the self-intersections of J- onto the (t,x)— plane are called shocks. 

In particular, we distinguish the genuine big front which is by definition the re- 
striction of the big front J 7 to {x G dom(dv)}, that is to the points such that v is 
differentiate at Xq, with corresponding genuine branches and genuine shocks. 

Proposition 4.3. Suppose ~F(v) is a typical big front for the (HJ) equation, 

(1) If the genuine shocks strictly verify the entropy condition, then within suffi- 
ciently small time, the minmax solution preserves the shocks. As a consequence, 
it coincides with the viscosity solution. 
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(2) On the contrary, if some genuine shock violates the entropy condition, by this 
we include rarefactions then the minmax solution admits new shocks other 
than the genuine shocks 1 

Proof. 1) Let x be a singularity of v, x(t) the genuine shock generated from x, and 
M ± (t,x) the evolution of the two genuine branches whose intersection gives x{t)- Note 
pf := dxU^^t, We shall suppose that p$ < p$, the other case being similar. 

There is e > 0, such that there exists a neighborhood O e (0, e) x R of x{t) in the big 
front J 7 that contains exactly the two branches u as well as J~\ x \, where: 

(4.8) = { z {p) = (x + tH'(p),v(x) + tipH'ip) - H{p))) , p e H,P Q ]}- 

We claim that for e small enough, does not interfere with the selection of the minmax 
solution i?Qt>(x) for (t,x) G O. For this purpose, it is enough to show that there is no 
point in J 7 ?-, lying below the genuine shock (x(t), u(t, x(t)))> where 

X (t) = xf + tH'ipf), 
u(t, X (t)) = u^x^+tipfH'ip^-Hipf)). 

Suppose t (->• at G [pq ,Pq] is a C 1 function of t such that 

z{a t ) = {x + tH'(a t ), v(x) + t{a t H'{a t ) - H(a t ))) eJ*n{x = x(t)}. 

Define 

A(t) = (vT(0,Xt) +t{p-H\p-) -Hip'))) - (u(0,x)+t(a t H'(a t )-H(a t ))). 

Then A(0) = 0. Note that d x u~(0,x~j~) = d x u~(t,x(t)) because they lie in the same 
characteristics, 

A{t) = p- t x- + p^H'ip-) - Hip-) - a t H'(a t ) + H(a t ) + t(p; H" \p~)p- - a t H"(a t )a t ) 
= ptixit) - H\p-) - tH"{p-)p-) + PtH\p-) - H(p;) - a t H\a t ) + H{at) 
+tip^H"ipl)fi + a t (H'ia t ) - x(t))) 
H(a t )-Hfo). 



(p t - a t)(x(t) 



<*t-Pt 



If x{t) strictly verifies the entropy condition, then we have a := lim t> o j t_>.Q a t < p . 
Indeed, if «o = p , then 

X(0) := Jim o x(t) = jm^'K) + tH\a t )a t = H\a Q ) = H\p ), 

which means that the shock x does not strictly verify the entropy condition. Therefore, 
for sufficiently small time t, we can assume that a t < p~t , then by the entropy condition 
(4.7), we get A{t) < 0, hence A(t) < 0. Now the only way to select a continuous section 



tr rhe rarefaction means there is separation of genuine branches. 
*By this, we may include the disappearance of the genuine shocks. 
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in J-^lo is choosing the two branches u with the genuine shock x{t)- Since the graph 
of the minmax solution {(x, Rqv(x))} is a continuous section in J 7 *, we get 



Rqv{x) 



u + (t,x), (t,x)eO + ux 

u~(t,x), (t,x) G 0~ U x 

by the local property of the minmax. 

2) Suppose that a genuine shock violates the entropy condition: if its two branches 
separate, there must be a piece from J^-y to compensate in order to get a continuous 
section; if they intersect, the minmax solution can not preserve the genuine shock, 
otherwise it satisfies the semi-group property, from which, followed by Proposition 3.1, 
that it must be the viscosity solution, thus a contradiction. □ 

Example 4.4. If H is convex (resp. concave), then in the wave fronts, the branches 
from the singularites of the initial function always lie above (resp. below) the the 
genuine branches, hence the min (resp. max) solution is the viscosity solution. See 
Figure 2. 




Figure 2. 



Example 4.5. If if is non convex, but the genuine shock does not violate the entropy 
condition, the configuration of the wave front are depicted in Figure 3. 




Figure 3. 



Before a critical time t c where the genuine branches violate the entropy condition, the 
graph of the viscosity solution, composed of the genuine branches, is always contained 
in each wave front ^(v), hence in the big front J~(v). After that, the behavior of 
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the viscosity solution becomes more complicated since the branches generated by the 
singularities of the initial function begin to interfere. As a consequence, the viscosity 
solution may not be contained in the big front. 

Instead of considering the big front U t > tc {t} x ^(v), let us consider 

T = U tem {t}xT t (R^v) 

for e > small. 

Let v c (x) := R l c v(x), suppose that RqV c (x) has a C 1 shock x{t) in (t,x) G O e : = 
(0, e) x O and that x(t) violates the entropy condition. Given a subdivision {( n } 
of (0, e), where ( n = {0 = t < < ■■■ < t^ = e}, let Xn(t) be a continuous, 
piecewise C 1 shock of the related iterated n-step minmax u n (t,x) := R,Q n v c (x). Note 
Pt := 9 x u n (t, Xn{t)±). By Proposition 4.3 and Lemma 3.10, for \( n \ small enough, 
there is violation of the entropy condition: for t G (t^ ^t+i]i Xnif) is constructed 

by the intersection of two branches, one of which comes from T k (n) (u n (tj^ , •)). 
Assuming that 

Xn(t) = Xn(tP) + (t~ tP)H'(p-), t G 

where p t lies between p\ n) and p~ {n) for t G (4™^am-i]' no ^ e that Xn(^) satisfies the 

^k ^k 

Rankine-Hugoniot condition, 



Xn(4" } ) = \ =-±- = H'{p- n) ), < k < n - 1. 



P. (n) P.(n) "* 

Thus the limit RqV c (x) of the sequence of iterated minmax {R l Q n v c (x)} n has a C 1 
shock satisfying 

(«) i(t) = " = tffe-), 0< ( <e. 

Pt - Pt 

By Theorem 3.13, the limit -Rq^^) is the viscosity solution, hence (4.9) describes 
the propagation of singularities of the viscosity solution. 

We see that those characteristics which are not contained in the original big front 
U t> t c {t} x ^(v) are tangent to the shock x(£), explained by (4. 9). This kind of shock is 
called a contact shock. The contact shock is a typical phenomenon due to rarefaction. 
With the limiting process of iterated minmax, we can see how arise the "exterior" 
characteristics: they are indeed emerging from the singularities of the minmax iterated 
step by step. 

In their paper [IK96], S.Izumiya and G.T.Kossioris have studied the bifurcation 
of shock waves for the {HJ) equations in the neighborhood of the critical time t c 
for generic Hamiltonians and initial functions. By using the method of Legendrian 
unfoldings, they classified the local pictures of the bifurcations of the viscosity solution 
into eight generic cases. 
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Here we will not repeat the description of the generic cases but take several concrete 
examples to explain the limit process of iterated minmax and hence the singularities of 
viscosity solutions after the critical time where there are violations of entropy condition 
for the genuine branches. 

Let x be a singularity of v , we first figure out how the wave front looks like by 
indicating the relative position of J^j-y = J-\{ Xo =x} in J 7 = Indeed, we get from 

(4.8) some key rules to characterize F\%y 

(1) The cusps in .FL, correspond to the inflection points of H; and the convexity 
of the branches of Fh* coincide with that of H; 

(2) The tangency of J-}g, at z{p) is p; 

(3) Fjz} and the two genuine branches are joined in a C l smooth manner at the 
extremities z(pq) with tangency p$ . 



Example 4.6 (Rarefaction). Consider 





Look at a neighborhood of the singularity x = of v. 





V 



H 



For t > s > small, the wave fronts are depicted as follows 
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In the wave front ^(v), the two branches in blue are genuine branches, and the 
curve in red is J-| |(t>). There is rarefaction, that is, separation of genuine branches. 
The 1-step minmax R\ ) v{x) attains the minimum in the wave front, and the shock 
Xi(t) of Rqv(x) is given by the intersection of the right genuine branch and a branch of 
J-| |(f). Revealed by the difference of ip f (dv) and ip t ~ s (dRQv), the shock Xiif) violates 
the entropy condition for t small, hence, as depicted in the wave front J zt ~ s {R s Q v), there 
is rarefaction for the genuine branches of RqV and the 2-step minmax R l s o R^v{x) has 
a shock X2{t) which, for t > s, is given by the intersection of the right genuine branch 
of Rqv(x) and a section of By the limiting process of iterated minmax, 

we then get a contact shock as described in the left of Figure 13. 

Example 4.7. Consider 

\x[x + 1), x < (J 
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Look at a neighborhood of the singularity x = of v: 





Figure 11. <p\dv) Figure 12. ^- s (dR s Q v) 
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As depicted in the wave front •7 r *(i>), the 1-step minmax R^v (x) has two shocks Xi(t) 
and Xi{t)- Both violate the entropy condition, hence the 2-step minmax Rl o Rqv(x) 
has new shocks xiW and X2W for t > s. Applying the same argument as before, 
by the limiting process of iterated minmax, we get that the viscosity solution has two 
contact shocks. See the right of Figure 13. 




Figure 13. 
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